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Abstract 

We show that the fermion determinant for 2-D Wilson lattice fer- 
mions coupled to an external scalar field is equivalent to self avoiding 
loops interacting with the external field. In an application of the re- 
sulting formula we integrate the scalar field with a Gaussian action to 
generate the A^-component Gross-Neveu model. The loop representa- 
tion for this model is discussed. 
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1 Introduction 



The fermion determinant is a highly non-local object as can e.g. be seen from 
the hopping expansion for lattice regularized fermions (see Q for a basic in- 
troduction). The hopping expansion expresses the fermion determinant as 
the exponential of a sum over all possible closed loops and the external field 
variables along a loop are collected as factors for this loop. The exponen- 
tial function can then be expanded and the result is a representation of the 
fermion determinant in terms of loops. The loop-sum in the exponent, how- 
ever, contains loops of arbitrary length, which can also iterate parts or all of 
their contour arbitrarily often. Thus, the exponent contains arbitrarily high 
powers of the external fields. On the other hand we know that, at least on a 
finite lattice, the fermion determinant is a finite polynomial in the external 
fields, and thus infinitely many contributions present in the exponent have 
to cancel each other when expanding the exponential. The result has to be 
a relatively simple loop representation for the fermion determinant. 

For the case of staggered fermions several papers can be found in the lit- 
erature ||2|-(5| where polymer representations for the fermions are obtained. 
For the case of Wilson fermions relatively few is known due to the more in- 
volved spinor structure of the fermions. Here we concentrate on 2-D Wilson 
fermions. 

An instance where the above discussed cancellation of contributions was 
brought under control for a model with Wilson fermions is Salmhofer's map- 
ping of the strongly coupled lattice Schwinger model to a self avoiding loop 
model p. In a first step the gauge fields at strong (=infinite) coupling were 
integrated out and the remaining Grassmann integral was then represented 
as a sum over loops. Subsequently Scharnhorst studied the two-flavor lat- 
tice Schwinger model with this method Q and extended the techniques to 
flnd a two-color loop model for the 2-D lattice Thirring model Q. In all 
of these cases, however the external field was either integrated out in the 
strong coupling limit ||6|, 0] or wasn't present at all |^]. On the other hand, 
our arguments given in the first paragraph show that the cancellation of 
higher winding loops is a universal phenomenon and that it should be pos- 
sible to find a simple loop representation also for the fermion determinant 
in an external field. 

In this article we present a simple formula for the fermion determinant 
of 2-D Wilson lattice fermions in a scalar background field (Eq. (|l2|)). The 
proof is based on a result for the hopping expansion for a generalized 8- 
vertex model where the vertices are coupled to a background field . With 
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a proper choice of the vertex weights the square of the partition function for 
this model represents the 2-D fermion determinant in a scalar background 
field. The resulting expression reduces the hopping expansion to a finite 
sum (on a finite lattice) of loops. In addition it is possible to explicitly 
integrate out the external field. By doing so with a Gaussian action for 
the scalar we generate a loop representation of the 2-D Gross-Neveu model. 
Such loop representations for lattice field theories allow for a considerably 
more accurate numerical treatment as has e.g. been demonstrated for the 



strongly coupled Schwinger model |10 |. 



2 Setting and hopping expansion 

The basic idea for the proof of our loop representation is to identify the 
hopping expansion of the Wilson fermion determinant with the hopping ex- 
pansion of a generalized 8- vertex model Q. Here we briefly rederive the 
hopping expansion for Wilson fermions in a form suitable for comparison 
with the generalized 8-vertex model (see e.g. [|| for a more detailed discus- 
sion of the hopping expansion). 

We study a lattice model of 2-D fermions which in the continuum corre- 
sponds to the action 

5 = /d^.«.) [.A (1) 

where 9{x) is a scalar external field. The lattice fermion determinant is 
expressed as a path integral 

detM[e] = f D^pDlp expl- ^ lp{x)M[e]{x,y)ij{y)\ , (2) 

V x,yeA J 

where the kernel for the lattice fermion action (Wilson fermions) which 
regularizes the continuum action (||) is given by 

±2 

M[e]{x,y) = [2 + m + 9{x)'\ 6^^y - ^ r^<5x+/.,j/ , 

M=±i 

and we defined 
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Here cri, (T2 are Pauli matrices. The sum in the exponent of (|2|) runs over the 
whole lattice A, which for simplicity we assume to be a finite rectangular 
piece of H? (the generalization to e.g. a torus is straightforward). The 
boundary conditions are open, i.e. hopping terms that would lead to the 
outside of our lattice are omitted. We define 

h{x) = 2 + m + B{x) , 

and assume that Q{x) is such that h{x) 7^ for all lattice points x. This is 
a purely technical assumption due to the particular techniques we use for 
computing the determinant. The final result will be a finite polynomial in 
the Q{x) and the above restriction is irrelevant then. We now can write 

detM[0] = ^h(xf^&i(\-R{0^ = \{h{xfeiivl-f^-TiR[e\A . 

xgA xeA \ n=l ^ I 

(3) 

In the last step the hopping expansion was performed, i.e. the determinant 
was expressed using the well known trace-logarithm formula and the loga- 
rithm was expanded in a power series. The hopping matrix R\0\ is defined 
as 

R\e\{x,y) = T,^^5.+,,y. (4) 

The series in the exponent of (|3|) converges for < 1, which can be 

enforced by choosing large enough i.e. suitable 9{x). Again this is 

only a technical restriction and can be abandoned in the final result. Due 
to the Kronecker delta in (^, the contributions to Tri?[0]" are supported 
on closed loops on the lattice, and since closed loops are of even length, the 
contributions for odd n vanish. For even n = 2k we obtain 

Tri?[^P = E E n Trnr,. (5) 

(2k) 

Here Cx is the set of all closed, connected loops of length 2k and base 
point X. By P{1) we denote the set of all sites visited by the loop I. Note 
that a factor l/h{x) is produced whenever I runs through x which can be 
arbitrary often for long enough loops. The last term in (^) is the trace of 
the ordered product of the hopping generators T^ as they appear along the 
loop /. We remark, that r-|-^r=|i^ = 0, which implies that whenever a loop 
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turns around at a site and runs back along its last link this contribution 
vanishes. Thus all these back-tracking loops can be excluded from 

Evaluating the trace over the matrices for a given loop is the remain- 
ing problem in the hopping expansion. It first has been solved in by 
realizing that the Pauli matrices give rise to a representation of discrete ro- 
tations on the lattice. Alternatively one can decompose the loop using four 
basic steps and compute the trace in an inductive procedure along the lines 
of 1 12]. The result is 

TrRr. = -(-lr«(-^)^'^ (6) 

By s{l) we denote the number of self-intersections of the loop / and c{l) gives 
its number of corners. The result is independent of the orientation of the 
loop. Inserting (^) and (^) in @ we obtain 



detM[^] = nM-)^exp E^E E i-'y'K\t' n jt) 

(7) 

Finally we further simplify this expression by removing the explicit summa- 
tion over the base points x. A loop of length 2k without complete iteration 
of its contour allows for 2k different choices of a base point thus cancelling 
the factor 1/2/;; in (0). A loop which iterates its whole contour I {I) > 1 times 
allows only for 2k/ 1(1) different base points and a factor !//(/) remains. The 
final expression is 

where C is the set of all closed, connected, non back-tracking loops of arbi- 
trary length. Each loop is included in C with only one of its two possible 
orientations and we collect an overall factor of 2 in the exponent. 



3 Identification of the corresponding generalized 
8- vertex model 

As already outlined, the next step is to compare the result (^) to the hopping 
expansion for a generalized 8- vertex model studied in detail in |^]. In this 
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generalized model, the vertices are coupled to a locally varying external 
field f{x). Before we proceed let's first discuss this generalized model and 
its relation to the standard 8-vertex model. 



The standard 8-vertex model |13, 14, ^ can be viewed as a model of 8 



quadratic tiles (vertices) and each of them is assigned a weight Wi {i = 1, ...8) 
(compare Fig. |l|). A tiling of our lattice A (the same rectangular piece of 
TL^ as before) is a covering of A with the tiles such that on each site of A we 
place one of our tiles with the centers of the tiles sitting on the sites. The 
set T of admissible tilings is given by those arrangements of tiles where the 
black lines on the tiles never have an open end. The partition function of 
the standard 8-vertex model is the sum over all admissible tilings t £ T and 
the Boltzmann weight for a particular tiling t is given by the product of the 
weights Wi for all tiles used in this tiling t. 
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Figure 1: The eight vertices (tiles) and their weights Wi. 

In our generalization of the model we now couple the vertices to an 
external field (p{x) located on the sites x of A. The partition function of the 
generalized model is given by 



n 



ter i=i 



ni{t) 



n ^(^) 

x£P{t) 



(9) 



Here P{t) denotes the set of all sites occupied by the tiling t. When a site x 
is occupied by tile Nr. 2, this site is counted twice giving a factor tp{x)'^ . In 
case X is occupied by tile Nr. 1, x ^ and the factor is 1. For all other 
tiles X is counted once and the factor is ^{x). 

The generalized model (P) also allows for a hopping expansion which is 
derived in The central step is to rewrite the partition function as an 
integral over Grassmann variables along the lines of |16, 17, O]. The action 
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is a quadratic form in the Grassmann variables and the partition function 
gives rise to a Pfaffian. The Pfaffian however can be expanded similar to the 
expansion in Section 2 and it is furthermore possible to explicitly evaluate 
the traces over the hopping matrix. The resulting expression reads g] 

zM = (-i)i^i-r' exp(E^^(^)'"n-r^'^ n ^i-) 

\1&C i=3 x&P{l) 

(10) 

As for the Wilson fermions the sum runs over the set C of all closed, non 
back-tracking loops of arbitrary length. Each loop I is included in C with 
only one of its two possible orientations. By |Z| we denote the length of the 
loop, /(/) is the number of iterations of its complete contour and s{l) is the 
number of self-intersections. P{1) again denotes the set of sites visited by 
the loop / with x being included in P{1) whenever the loop I runs through 
X. By |A| we denote the size of the lattice. In the following we will work on 
lattices with an even number of sites and the overall sign factor is 1. 

The exponents ni{l) in ( p!0| ) give the numbers for the abundance of the 
line elements as they are depicted in Fig. ||. E.g. when the loop / changes 
from heading east to heading north at a site it picks up a factor of and 
similarly for the other tiles w^,W4^,w%,wt,w%. Note that the loops / G C 
occur as an ordered set of instructions for the directions the loop takes as it 
hops from one site to the next. The element corresponding to the tile 2 with 
weight W2 (compare Fig. |^) is not needed to describe the loop /. Equation 
(|lO|) does not contain W2 explicitly at all. The weight W2 is related to the 
other weights through the free fermion condition [p!^ ] 

U1UJ2 + UJ^UJi = + ■ (11) 

The free fermion condition is a sufficient condition for finding an explicit 
solution to the standard 8-vertex model without external fields. The above 
mentioned Grassmann representation automatically enforces the free fermion 
condition as is discussed in more detail in [^. 

By comparing ^ with ( p!o[ ) we find that the sums over the loops in the 
exponent become identical when setting 

1 , , 1 

Wl = = Wi = 1 , = Wq = Wt = = —7= , ^\x) 



Using the free fermion condition (11) we find W2 = This implies that when 
representing Wilson fermions using Eq. (0), the sum over the tilings can be 
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replaced by a sum over the set Lsa of closed, self-avoiding loops, i.e. loops 
that are not allowed to self intersect or touch each other. The loops can 
have several disconnected but closed pieces and each piece is included in Lsa 
with only one of its two possible orientations. 

It is important to note, that the exponents in (^) and ( p!o[ ) differ by an 
overall factor of 2. This is due to the fact that the action for the Wilson 
fermions is a bilinear form giving rise to a determinant while the Grassmann 
action for (^) is a quadratic form giving rise to a Pfafhan when integrating 
out the Grassmann variables. Since the kernel of the latter action is anti- 
symmetric the Pfaffian, however, is given by the square root of a determinant 
causing the difference by a factor of 2 in the exponent. 

Putting things together we obtain the final formula 

..M„.np....,.„^ (_£(-L)-n_^_L_ 

(12) 

where as above c(^) is the number of corners of I. The loops in Lsa ai"e 
self-avoiding and thus for a given loop configuration / each site of the lattice 
is occupied only once. Thus in the sum each inverse field [2 + m + 
can only occur linearly. When taking the square of this sum we can only 
produce terms which are at most quadratic in the inverse field. The overall 
factor, however, still can cancel the quadratic terms and the final result is a 
finite polynomial in the fields Q{x). The above imposed technical restrictions 
on the range of Q{x) can now be lifted. 

Several remarks on the result ([T^ are in order: When setting all external 
fields to zero we find that free Wilson fermions are equivalent to the square 
of the partition function of the self avoiding loop model ||T9| , pO| with bending 
rigidity l/\/2 and bond weight [2 + ?7i]^^. Thus for a trivial background field 
we reproduce the result for free fermions obtained by Scharnhorst with 
a different method. Eq. (|l2|) is a direct generalization of the trivial case to 
fermions in a scalar background field. 



From an algebraic point of view Eq. (12) is exactly the expression we had 
in mind when discussing the general algebraic structure of the determinant. 
In two dimensions the fermion determinant has to be a polynomial which 
can at most be quadratic in the external field. The feature that the terms 
in this polynomial are organized according to closed loops is inherited from 
the hopping expansion. Actually, for single contributions to (O) it is even 
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possible to trace their emergence from an expansion of the exponential in 
(P). When doing so, one finds that the intersection factor (— l)'^^') provides 
the mechanism which ensures the cancellation of loops with multiply oc- 
cupied links. Eq. (|lO| ) establishes that this cancellation mechanism is also 
independent of the corner weights. 

Finally, it is obvious from (jl^) that in the loop representation it is 
straightforward to integrate out the scalar fields. This can be done with 
different actions for 0. In the next section we discuss the case of a simple 
Gaussian which will produce the Gross-Neveu model |21]. We remark, that a 
loop representation for the Gross-Neveu model with staggered fermions has 
been analyzed in O and a numerical study of the model with conventional 



methods (introduction of the auxialiary field) is given in [ 22 



4 Application of the result to the Gross-Neveu 
model 

In order to give an application of our Eq. (^) we now integrate the scalar 
field with a Gaussian measure 



dii[e] 



n 



de{x) 



f2^ 



exp ( -^^(^ 



When integrating (detM[0])^ with this measure we generate a lattice ver- 
sion of the A^-component Gross-Neveu model pTl with (continuum) action 



S 



d^x{^lj{x) 7^3^ — m il){x) — il){x)ip{x) 



(13) 



where ijj^ijj now have N flavor components. Using (|T^) the lattice partition 
function for the Gross-Neveu model reads 



2N 



dfi[9] n [2+m+e{x)] 



2N 



^ (^)'" n 



h,---,hN 



E/ 1 



1 ^c{«i)+...+c(/2iv) 



yeP{l) 



2+m+9{y) 



l[G[2N-0,{h,...,lN)]. (14) 
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The function Ox{h, I2N) counts how many of the independent loops h, I2N 
occupy the site x for a given loop configuration. Each loop can either leave 
the site empty or occupy it once (the loops are self-avoiding) and thus Ox- 
has values between and 2N. The function G[J], J = 0, 1, 2N is given 
by 



G[J] = J^^~^^'^i^+^+(^V = E'(o^J(2^-l)"/(2 + 




Equation (|T^ establishes that the A^-component Gross- Neveu model is equiv- 
alent to a model of 2N independent self-avoiding loops. The partition func- 
tion is a sum over all loop configurations with a weight corresponding to the 
total number of corners times a simple function of the occupation number 
for each site. 

The partition function can easily be reformulated as a 8^^-vertex model. 
The new tiles are obtained by using 2N different colors and with each color 
one of the line patterns of Fig. || is drawn onto our new tile, giving a total 
of 8^^ different tiles. The weight for the tiles is given by a product of the 
weights Wi for each color (1 for 'Wi,W3,W4, for W2 and l/\/2 for w^, ...,ws) 
times the occupation factor G[J]. 

Let's discuss the case of = 1 in more detail. Here the tiles show lines 
in two colors, say red and blue. The weight factors are the products of the 
Wi for the blue and red lines multiplied by 1 when the tile is empty, by a 
factor of 2 + m when there is only one color and a factor of (2 + m)^ + g 
when both colors are used on the tile. When setting g = the weights can 
be factorized into two terms corresponding to the two colors and we recover 
the result for the free case already discussed above. 

A second choice of parameters also leads to a particularly simple model. 
When setting 2 + m = 0, we find that all tiles which show only one color 
vanish. Thus for a non-vanishing contribution every site has to either be 
empty or is visited by both, a blue and a red loop. For our simple rectangular 
lattice this implies, that the red and blue loop configurations have to sit on 
top of each other and integrating out the scalar field at 2 + m = has 
produced a Dirac delta on the space of loops. The resulting model is again 
a self avoiding loop model, now with bending rigidity 1/2 and bond weight 
g. This is the same model which was shown |^] to be equivalent to the 
strong coupling Schwinger model. It has been analyzed with Monte-Carlo 
methods in [11^] where the existence of a phase transition at gc = 0.5792 
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was established. For a study of this model, using computer algebra on small 
lattices see [23|. In the form of the above discussed two-color (64- vertex) 
model, the = 1 Gross-Neveu model is now relatively simple to analyze 
numerically in the whole g, m-plane. In particular, the representation as a 
vertex model may allow for the use of efficient cluster algorithms |24]. 



5 Summary and discussion 

The motivation for this article was to find simple representations for the Wil- 
son lattice fermion determinant in an external field. Counting the powers 
of the external field in the determinant one finds that a large set of contri- 
butions appearing in the standard hopping expansion has to cancel. The 
determinant can only be a finite polynomial in the external variables. Find- 
ing, however, a reasonably simple and useful expression for the determinant 
is a hard problem. 

In this article we succeeded in finding such a simple representation for the 
case of 2-D Wilson fermions in a scalar background field. The determinant 
can be written as the product of two self-avoiding loop models coupled to 
the external field. To prove this result, Grassmann techniques and hopping 
expansion for a generalized 8-vertex model were used. In the obtained loop 
representation it is straightforward to integrate out the scalar fields, and 
the application to the case of the Gross-Neveu model was discussed in more 
detail. 

When one instance of a simple expression for a Wilson type fermion 
determinant in an external field can be found it is natural to ask whether 
other and more realistic models also allow for a simple representation of the 
fermion determinant. 

The next candidate are 2-D fermions interacting with an abelian vector 
field. The vector field is coupled to the fermions through its gauge trans- 
porters which are supported on links instead of sites. This does not pose 
a fundamental problem since the generalized 8-vertex model Q can also 
be formulated for external field living on links. All terms of the hopping 
expansion can still be evaluated explicitly. A certain difficulty is, however, 
given by the fact that the link variables have to be complex conjugated when 
hopping in backward direction. Thus the loops in the hopping expansion 
give contributions complex conjugated to each other when the orientation is 
reversed. We believe, however, that when properly organizing the loops this 
problem can be overcome. A hint in this direction is Scharnhorst's proof 
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for the existence of a vertex model for the 2-D Thirring model Q . If a loop 
representation for the determinant in an external vector field exists then it 
should be possible to obtain Scharnhorst's result by integrating out the vec- 
tor field, similar to the loop representation of the Gross-Neveu model which 
was obtained by integrating out the scalar field. The case of non-abelian 
gauge fields poses the additional difficulty, that due to the traces over the 
gauge field matrices the product of two loops cannot be written as a new 
loop. Nevertheless, the power counting argument discussed in the introduc- 
tion still suggests, that representations simpler than the standard hopping 
expansion should exist. Similar in 4 dimensions where at least the problem 



of computing the trace of the 4-D 7-matrices is solved [11 
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